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ABSTRACT 

We report the first detection of X-ray emission from a brown dwarf in the Pleiades, the M7- 
type Roque 14, obtained using the EPIC detectors on XMM-Newton. This is the first X-ray 
detection of a brown dwarf intermediate in age between w 12 and ~ 320 Myr The emission 
appears persistent, although we cannot rule out flare-like behaviour with a decay time-scale 
> 4 ks. The time-averaged X-ray luminosity of Lx ~ 3.3 ± 0.8 x 10^^ erg s^^, and its ratios 
with the bolometric (-Lx/iboi ~ lO"'^"^) and Ha (Lx/Lna ~ 4.0) luminosities suggest 
magnetic activity similar to that of active main-sequence M dwarfs, such as the M7 old- 
disc star VB 8, though the suspected binary nature of Roque 14 merits further attention. No 
emission is detected from four proposed later-type Pleiades brown dwarfs, with upper limits 
to Lx in the range 2.1-3.8 xlO^'' ergs^^ and to log(Lx/-Z^boi) in the range —3.10 to —2.91. 
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1 INTRODUCTION 

Magnetic activity, generating chromospheric Ha and coronal X-ray 
and radio emissions, is a ubiquitous feature of main-sequence (MS) 
late-type stars (spectral types ~F5-M7). Studies of these diagnos- 
tic emissions have found consistency in the character of magnetic 
activity throughout this range, despite the expected change in dy- 
namo mechanism demanded by the absence of a radiative interior 
in stars of spectral types ~ M3 and later However, recent studies 
suggest the magnetic activity of 'ultracool' objects, with spectral 
types ~ M8 and later, is quite different. 

The observed persistent ('non-flaring') levels of X-ray 
(Lx/Lhoi) and Hq (Lua / Lhoi) emission from MS late-type stars 
increase with decreasing Rossby number, Ro = P/tc, where P 
is the rotation period and re is the convective turnover time, un- 
til reaching respective 'saturation' plateaus of Lx/iboi ~ 10~^ 
and Luc/Lboi ~ 10"^'^ (e.g. Delfosse et al. 1998 for M dwarfs). 
The fraction of field stars showing persistent chromospheric emis- 
sion levels close to saturation increases toward later spectral types, 
peaking around M6-7 (Gizis et al. 2000). 

However, around spectral type M9 persistent Ha emission lev- 
els begin to plummet dramatically (Gizis et al. 2000). Among L- 
type dwarfs no rotation-activity connection is found: Lua/Lhoi 
continues to fall steeply toward later spectral types despite most L- 
dwarfs being fast rotators (Mohanty & Basri 2003). A proposed ex- 
planation is that magnetic fields diffuse with increasing efficiency 
in the increasingly neutral atmospheres of cooler dwarfs (Meyer 
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& Meyer-Hofmeister 1999; Mohanty et al. 2002), overwhelming 
the importance of a rotation-driven dynamo efficiency in chromo- 
spheric heating. Magnetic activity is still observed, however, in the 
forms of Ha flaring on some L dwarfs (e.g. Liebert et al. 2003), 
and flaring and apparently-persistent radio emission from several 
ultracool dwarfs (Berger et al. 2001; Berger 2002). 

Detections of X-ray emission from ultracool field dwarfs are 
scarce. The M7 old-disc star VB 8 shows persistent and flaring 
X-ray emission levels of Lx/iboi ~ 10~^'^-10~^*, similar to 
those of active M dwarfs (Fleming et al. 1993; Schmitt, Fleming 
& Giampapa 1995; Fleming, Giampapa & Garza 2003). However, 
the persistent levels of X-ray emission from the M8 old-disc star 
VB 10 and the M9 -320 Myr-old brown dwarf LP 944-20 are at 
least an order of magnitude lower - Lx/Lhoi ~ 10~^" (Flem- 
ing et al. 2003) and Lx/Lboi < 10"^'^ (Rutledge et al. 2000), 
respectively - despite the latter being a fast rotator (v sin i = 30 
kms^^). Yet transient strong magnetic activity is evidenced by the 
flaring X-ray emission, with peak Lx/iboi ~ lO^^'^-lO^^'", that 
has been observed on both VB 10 (Fleming, Giampapa & Schmitt 
2000) and LP 944-20, and on the M9 field dwarfs LHS 2065 
(Schmitt & Liefke 2002) and IRXS J115928.5-524717 (Ham- 
baryan et al. 2004). Interestingly the temperature of the dominant 
X-ray-emitting plasma appears to be low, T ~ 10^'^ K, whether 
it is measured in the persistent (VB 10) or flaring (LP 944-20 and 
IRXS 1115928.5-524717) emission state. While such low temper- 
atures are typical for the persistent coronae of inactive stars - M 
dwarfs (Giampapa et al. 1996) and the Sun (Orlando, Peres & Reale 
2000) alike - the temperatures of flaring plasma are significantly 



higher, with T > 10' 
lando 2001). 



K (Gudel et al. 2004; Reale, Peres & Or- 
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2 Briggs & Pye 

Table 1. Observed and derived physical parameters of the studied brown dwaii's. Columns show (1-3) Name and alternative designations; (4-5) RA and Dec 
(J20()0; Pinfield et al. 2000); (6-8) J, H and K magnitudes (Pinfield et al. 2000, i-cr uncertainties are < 0.05 mag); (9) spectral type; (10) reference for 
spectroscopy: Rebolo et al. 1995 (R95), Zapatero Osorio et al. 1997b (Z97), 1999 (Z99), Martin et al. 1998 (M98); (11) effective temperature in K (assigned 
from spectral type after Mohanty & Basri 2003); (12-13) mass (Mq) and bolometric luminosity (10'"' ergs~^) calculated from T^g using models of Baraffe 
et al. (1998) for age 125 Myr; (14) binary mass ratio where binarity is indicated (Pinfield et al. 2003). 



Name 
(1) 



BPL 

(2) 



NPL 

(3) 



RA 

(4) 



Dec 

(5) 



J 

(6) 



H 

(7) 



K 

(8) 



SpT 
(9) 



Ref 

(10) 



Toll 
(11) 



(12) 



Lbal 

(13) 



9 
(14) 



Roque 14 


108 




03 46 42.90 


+24 24 50.4 


15.50 


14.89 


14.47 


M7.0 


Z97 


2700 


0.066 


3.7 


0.75-1 


Roque 12 


172 


36 


03 48 18.96 


+24 25 12.9 


15.91 


15.36 


15.10 


M7.5 


M98 


2625 


0.061 


3.1 




Teide 1 


137 


39 


03 47 17.83 


+24 22 31.5 


16.32 


15.62 


15.08 


M8.0 


R95 


2550 


0.056 


2.7 


0.5-0.9 


Roque 1 1 


132 


37 


03 47 12.02 


+24 28 31.5 


16.03 


15.64 


15.13 


M8.0 


Z99 


2550 


0.056 


27 




Roque 9 


100 




03 46 23.05 


+24 20 36.1 


16.28 


15.63 


15.22 


(M8.0) 




(2550) 


(0.056) 


(2.7) 





As very young substellar objects (f ^ 5 Myr) may have pho- 
tospheres as warm as MS M5-6 dwarfs, an individual brown 
dwarf may experience a transition from 'stellar-hke' to 'ultra- 
cool' magnetic activity as it cools. Brown dwarfs in star-forming 
regions are routinely observed to emit X-rays at high levels, 
ix/iboi > 10"^-*, arising from plasma at T> lO'^" K, similar 
to those of dMe stars and higher-mass young stars (e.g. Neuhauser 
& Comeron 1998; Imanishi, Tsujimoto & Koyama 2001; Mokler & 
Stelzer 2002; Preibisch & Zinnecker 2002; Feigelson et al. 2002). 
The ~ 12 Myr-old, low-mass brown dwarf TWA 5B, of ultracool 
spectral type M8.5-9, exhibits apparently persistent X-ray emis- 
sion at I/x/iboi ~ lO""**, like younger brown dwarfs, but with 
T « 10®-^ K (Tsuboi et al. 2003), like field ultracool dwarfs and 
LP 944-20. 

At around spectral type M8-9, between the ages of ~ 10 and 
~ 300 Myr for brown dwarfs, persistent X-ray emission levels ap- 
pear to fall by a factor ~ 100 and coronal temperatures appear 
constrained to T < lO'''^ K, even during flares and at high emission 
levels. The population of brown dwarfs in the Pleiades cluster, 135 
pc away (Pan, Shao & Kulkami 2004), of age ^ 125 Myr (Stauf- 
fer, Schultz & Fitzpatrick 1998) and spanning spectral types M6.5- 
early-L (Martin et al. 1998), is therefore crucial to understanding 
the evolution of substellar magnetic activity and the conflict of a 
rotationally-driven magnetic dynamo against atmospheric neutral- 
ity. ROSAT observations of brown dwarfs in the Pleiades detected 
no X-ray emission at the level of Lx/iboi ^ 10^^'^ (Neuhauser 
et al. 1999, and see Section 4.4). 

We present a deeper X-ray (0.3-4.5 keV) observation of five 
candidate brown dwarfs in the Pleiades (described in Section 2), 
using the the more sensitive XMM-Newton observatory (Section 3). 
We detect X-ray emission from the M7 Roque 14, investigate its 
temporal and spectral nature, and place upper limits on the X-ray 
emission levels of the undetected brown dwarfs: Teide 1, Roque 9, 
Roque 11 and Roque 12 (Section 4). We discuss the relative X- 
ray and Ha emissions of these objects in the context of magnetic 
activity on ultracool dwarfs and its evolution (Section 5) and close 
by summarising our findings (Section 6). 



2 SAMPLE OF BROWN DWARFS IN THE PLEIADES 

Five objects in our XMM-Newton field have been proposed as can- 
didate brown dwarf members of the Pleiades on the basis of optical 
and near-infrared (NIR) photometry in the IZJHK bands (Zapatero 
Osorio et al. 1997a; Pinfield et al. 2000; 2003). Observed and de- 
rived physical parameters are listed in TableQ All except Roque 9 
have published spectral types consistent with those expected of ~ 
125 Myr-old brown dwarfs. Further evidence for or against mem- 



bership of the Pleiades and hence substellar status is summarised 
below: 

Teide 1 is the on-axis target of the XMM-Newton observation. Its 
status as a brown dwarf member of the Pleiades is well-established 
on the basis of its proper motion (Rebolo, Zapatero Osorio & 
Martin 1995), the detection of Li in its spectrum and its radial ve- 
locity (Rebolo et al. 1996). It shows Ha emission with variable 
equivalent width, EWna = 3.5-8.6 A (Rebolo et al. 1995; 1996). 
It has been suspected to have a lower-mass companion from its po- 
sition on a JK colour-magnitude diagram (Pinfield et al. 2003). 
Roque 11 has an anomolous position on a JHK colour-colour 
diagram (Pinfield et al. 2003) but its radial velocity of —3.5 ± 7 
kms^^is consistent with those of other Pleiades members, and its 
Na I absorption is lower than that of field stars of the same spectral 
type, indicating lower gravity, and hence youth (Zapatero Osorio 
et al. 1997b). 

Roque 12 has a radial velocity consistent with Pleiades member- 
ship (Festin 1998), low Na I absorption and is a strong Ha emitter 
with EWna = 19.7 A (Martin et al. 1998). 

Roque 14 has low Na I absorption and strong Ha emission with 
EWna = 17.0 A (Zapatero Osorio et al. 1997b). It has been sus- 
pected to be a near-equal mass binary on the basis of its position 
in an IK colour-magnitude diagram (Pinfield et al. 2003) but no 
comparably bright companion with separation > 0.1 arcsec has 
been found (Martin et al. 2000). 

Roque 9 has no published spectral type but we estimate a spectral 
type of M8 as its NIR photometry is similar to those of Teide 1 and 
Roque 11. 

We shall use the term "brown dwarf" to refer to all five ob- 
jects, but note that the evidence in support of substellar status varies 
across the sample, and is weak for Roque 9. 



3 OBSERVATION AND DATA ANALYSIS 

The XMM-Newton observation, 0094780101, was centred on 
Teide 1 (J2000: a = 03''47'"18=0, 5 = +24°22'3l") and con- 
ducted on 2000 September I in orbit 134. The Thick optical block- 
ing filter was placed in front of all three EPIC detectors: the pn 
(Striider et al. 2001) was exposed for 40.6 ks and each MOS (Ml 
and M2; Turner et al. 2001) for 33.0 ks, beginning 7.5 ks later. 
The data were processed using the SCIENCE ANALYSIS SYSTEM 
(SAS) v5.4. 1^ and each EPIC eventlist was further filtered to ex- 
clude flagged 'bad' events, and uncalibrated event patterns (> 12 
for MOS ; > 4 for pn). Several short intervals affected by high back- 
ground were also excluded. We considered only events with PI in 
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Figure 1. (a) Adaptively smoothed EPIC 0.3^.5 keV image of the re- 
gion around Roque 14. The circle is centred on the NIR position, (b) Mo- 
saicked EPIC 0.3-1.4 keV image of a composite of the four undetected 
brown dwarfs. The dotted circle of 8-arcsec radius centred on the NIR po- 
sition marks the region used to estimate upper limits. The circled source 25 
arcsec to the west is not associated with the Pleiades. 



the range 300-4500 (nominally energies of 0.3-4.5 keV, and PI is 
reported in units of eV or keV from this point) to reduce back- 
ground contamination. 

We extracted an image with 4x4 arcsec square pixels from 
each detector and performed source detection in each image. The 
procedure, using tasks available in the SAS, is described in de- 
tail in Briggs & Pye, in preparation. In brief, potential sources 
were located using a wavelet detection package (EWAVELET), and 
masked out of the photon image while it was adaptively smoothed 
(using ASMOOTH) to generate a model of the background. The 
spatial variation of vignetting and quantum efficiency was mod- 
elled in an exposure map (produced using EEXPMAP). The images, 
background and exposure maps^ of the three EPIC instruments 
were also mosaicked to optimize the sensitivity of the analysis and 
EWAVELET was used to locate potential sources in the EPIC image. 
In each image, at the position of each EWAVELET source, a maxi- 
mum likelihood fitting of the position-dependent instrument point 
spread function (PSF) was applied (EMLDETECT) to parametrize 
each source and those with ML > 6 were retained. EMLDETECT 
additionally reconstructed a smooth model of the input image by 
adding PSFs to the background map with the location, normaliza- 
tion and extent of the parametrized sources. We repeated the source 
detection procedure on 20 images generated from this model image 
using Poisson statistics, and 20 images generated from the source- 
free background map, and found in both cases a mean number of 
false detections per image of 5 with ML > 6 and 2 with ML > 7. 

Approximately 130 X-ray sources were found, of which 34 
were associated with proposed members of the Pleiades. Match- 
ing of the EMLDETECT source positions with NIR positions of 
these members (Pinfield et al. 2000) revealed a boresight shift in 
the EPIC image of qx cos (5x — aNiR cos 5nir = 3.7 arcsec and 
5x — ^NiR = —0.8 arcsec, which is accounted for in the remainder 
of this work. The corrected positions of sources associated with the 
Pleiades were all within 6 arcsec of the NIR positions. The proba- 
bilities of a positional coincidence within 6 arcsec of a brown dwarf 
with a spurious detection or a source unrelated to the Pleiades are 
~ 10~^ and « 0.02, respectively. 

ML-fitting was performed in each image at the NIR position 
of each brown dwarf. The successful detection of Roque 14 is de- 
scribed further in Section 4. 1 , and investigations of the temporal 

^ The pn exposure map must be scaled by a factor to account for its higher 
sensitivity compared to Ml and M2. This factor is dependent on the source 
spectrum and discussed in Section 4.3. 
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Figure 2. (a) The arrival time and energy of each pn (black circle) and MOS 
(grey square) event detected in the Roque 14 source extraction region. In 
(b) the observation has been divided into four equal time intervals to bin all 
the events, (c) shows the expected contiibution of background events. The 
bold and dashed histograms are for events in the 0.3^.5 and 0.3-1.4 keV 
energy bands respectively. There is an apparent concentration of 0.3-1.4 
keV events in the second interval. 



and spectral nature of its emission are pursued in Sections 4.2 and 
4.3, respectively. In no other case was an ML value > 1 found. Up- 
per limits to the X-ray luminosity of these objects are calculated in 
Section 4.4. 



4 RESULTS 



4.1 An X-ray detection of Roque 14 



Roque 14 was returned as an X-ray source by EWAVELET in the 
M2, pn, and mosaicked-EPIC images. EMLDETECT determined ML 
values of 2.1, 5.9 and 7.0 (1.6, 3.0 and S.Scr) in these respective 
images, confirming the detection at ML > 6 (S.Ocr) only in the 
mosaicked-EPIC image'' (Fig. ^). The best-fitting X-ray source 
position in the EPIC image was 0.8 arcsec offset from the NIR po- 
sition. We estimated the I -a uncertainty in the relative positions 
as the sum in quadrature of the statistical fitting uncertainty, 0.73 
arcsec, an uncertainty in the EPIC absolute pointing of 1.0 arcsec, 
and an uncertainty in the optical position of 0.5 arcsec; thus the po- 
sitional offset of the source is ~ O.Gcr. The total number of EPIC 
source counts determined by EMLDETECT was 24.7 ± 5.9. We cal- 
culate the X-ray luminosity after consideration of the source spec- 
trum in Section 4.3. 

Further analysis of the X-ray emission froin Roque 14 has 
been conducted both including and excluding periods affected by 
high background, and no significant improvement has been found 
by excluding those periods. Therefore, to maximize the number of 
events available, we report the results of including the short periods 
of high background. 

We extracted source events from a circle centred on the best- 
fitting source position with radius 10.75 arcsec, which encloses 
only 54 per cent of the source counts (Ghizzardi 2002) but opti- 
mises the signalmoise ratio, and background events from the sur- 
rounding annulus, 8 times larger than the source extraction region. 

^ ML values output by EMLDETECT have been corrected as advised in 
XMMNews29. 
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10' 2x10' 3x10' 4x10' 1000 2000 3000 4000 
Time (s) PI (eV) 

Figure 3. Cumulative distribution functions of event anival times (a and 
c) and PI values (b and d). Bold lines are for observed background (grey 
in a and b) and source+background (black in c and d) data. Broken lines 
are for modelled source (a and b) and source+background (c and d) data. 
In (a) and (c) the dashed line is for a constant source count-rate, while dot- 
dashed and dotted lines assume all source counts arise from a flare with 
exponential decay on respective time-scales of 1000 and 5000 s. In (b) and 
(d) the dashed, dot-dashed, and dotted lines are for source temperatures of 
10^°, 10'^° and 10*° K respectively. A K-S test rules out the temporal 
model of a flare with decay time-scale 1000 s and the 10^" K spectral 
model at 90 per cent confidence. 



Ml and M2 eventlists were combined as these instruments have 
near-identical sensitivity, spectral response and exposure length. A 
total of 32 counts in all EPIC instruments (23 from pn and 9 from 
the two MOS cameras) was extracted from the source region and 
140 (99 pn; 41 MOS) counts from the background region. Thus we 
expect 17.5 (12.4 pn; 5.1 MOS) background counts in the source 
region. As the background region will contain ~ 46 per cent of 
the ~ 25 source counts, we expect a contribution of source counts 
(« 1.4) similar to the l-cr uncertainty (~ 1.2) in the estimated 
mean background counts. The probability of 32 counts or more ap- 
pearing in the source region as a result of a Poissonian fluctuation 
in a mean background of 17.5 counts is 0.0011 (3.2a), in support 
of the EMLDETECT detection. 

The events extracted from the source region were used to ex- 
amine the temporal and spectral behaviour of Roque 14. Events 
were extracted from a larger source-free region on the same CCD 
as Roque 14 to make a better model of the time and energy distri- 
bution of the background. 

4.2 Transient or persistent emission? 

Examination of the arrival times of the photons in the source 
eventlist suggests a concentration, particularly of lower-energy 
(0.3-1.4 keV) events, in the f» 10 ks interval 10-20 ks after the 
start of the pn exposure (Fig.|2j- 




0.4 0.1 0.2 0.3 0.4 

K-S value 



Figure 4. Confidence level of deviation from the test model as a function 
of K-S test value resulting from Monte Carlo simulations of (a) arrival time 
and (b) PI of 32 pn and MOS source and background events (see text for 
more details). The line-styles of black lines represent the same source mod- 
els as in Fig.|3] Unbroken lines show the standard calculated relation. When 
source and background distributions differ appreciably, the calculated con- 
fidence may significantly overestimate the true confidence. 



We investigated the statistical significance of this possible 
variability using the Kolmogorov-Smirnov (K-S) statistic. The K- 
S statistic is the maximum vertical difference between the cumula- 
tive distribution functions (CDFs) of an observed dataset and a test 
distribution (or two observed datasets). As our test model must ac- 
count for the significant number of background events in the source 
extraction region, which is subject to Poissonian fluctuations about 
its expected number, we have performed Monte Carlo simulations 
to assess the confidence level of deviation from the test model as a 
function of K-S value. 

Our test-model CDF was constructed from the arrival times of 
events in the background extraction region and arrival times of a 
number of source events (such that the ratio of source:background 
counts was as estimated from the observed data) chosen at uniform 
intervals in the CDF of a constant count-rate source. Sets of arrival 
times for pn and MOS events were constructed separately and then 
merged before calculating the test-model CDF for all EPIC data 
(Fig.H). 

Each trial dataset in our simulations was composed of the 
observed numbers of pn (23) and MOS (9) events. The number 
of background events was drawn from a Poisson distribution with 
mean value as estimated from the data. Arrival times for these back- 
ground events were drawn at random from the observed distribution 
for background events. Arrival times for the remaining events were 
drawn at random from the modelled distribution for source events. 
MOS and pn datasets were again constructed separately before be- 
ing combined. For each trial dataset the K-S statistic was calcu- 
lated. The CDF of these K-S values gives the confidence level of 
deviation from the test model as a function of measured K-S value 
(Fig.Eli). 

The simulated distribution of K-S values for a test model of a 
constant source count-rate, based on 10000 trials, is practically in- 
distinguishable from the standard calculated distribution. The K-S 
value of the observed dataset for this test model is 0.110 (Fig|3j:), 
which corresponds to just 20 per cent confidence that the source 
count rate is not constant. Following the same procedure we have 
additionally tested models in which the source emission is purely 
transient (Fig.|3ji), with the profile of fast (here, instantaneous) rise 
and exponential decay typical of flares on late-type stars, begin- 
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X-ray emission from a brown dwarf in the Pleiades 5 

Table 2. X-ray and Ha properties of the studied brown dwarfs. Columns show: (2) total MOS-equivalent effective exposure time in s; (3) number of EPIC 
source counts and 1-cr uncertainty, or 95 per cent upper limit; (4) MOS-equivalent count rate (10^'* s~^) and 1-cr uncertainty, or 95 per cent upper limit; 
(5) X-ray luminosity in lO^'^ ergs^^ in the 0.3^.5 keV band and l-cr uncertainty, or 95 per cent upper limit; (6) logarithmic ratio of X-ray and bolometric 
luminosities; (7) equivalent width of Ho (A); (8) logarithmic ratio of Ho and bolometric luminosities (calculated as in Mohanty & Basri 2003); (9) ratio of 
X-ray and Ho luminosities. 



Name 


T„,p 


JVx 


Cx 


ix 


l°g(ix/-£'boi) 


EW-a^ 


log{-C'Hc./ibol) 


ix/-C-Hc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Roque 14 


117590 


24.7 ± 5.9 


2.1 ± 0.5 


3.3 ± 0.8 


--3.05 


17.0 


-3.66 


4.0 


Roque 12 


26300 


< 6.5 


< 2.48 


< 3.8 


< -2.91 


19.7 


-3.69 


< 6 


Teidc 1 


194250 


< 26.9 


< 1.39 


< 2.1 


< -3.10 


3.5,8.6 


-4.53.-4.14 


< 27,< 11 


Roque 1 1 


145080 


< 21.1 


< 1.46 


< 2.3 


< -3.08 


5.8 


-4.31 


< 17 


Roque 9 


73640 


< 11.6 


< 1.57 


< 2.5 


< -3.04 









ning 10 ks after the start of the pn observation. For short decay 
time-scales, when the source and background models are starkly 
different (Fig.|3^), Fig.|4^ shows that the standard calculation sig- 
nificantly overestimates the simulated confidence level of deviation 
from the test-model. The K-S test rules out flare models with de- 
cay time-scale of 3 ks or less at the 90 per cent confidence level, 
but permits those with time-scale 4 ks or longer (Fig.|3j;). 

We therefore cannot exclude the possibility that the observed 
emission from Roque 14 was due to a flare-like outburst similar to 
those with decay time-scales of ~ 5 ks seen from the M9 ultra- 
cool dwarfs LP 944-20 (Rutledge et al. 2000), LHS 2065 (Schmitt 
& Liefke 2002), and RXS Jl 15928.5-524717 (Hambaryan et al. 
2004). 



4.3 Temperature of X-ray emitting plasma 

Fig. |2| suggests the excess of counts was chiefly in the low- 
energy (0.3-1.4 keV) band, indicative of plasma with temperature 
T < lO'^'^^ K typical in coronae on late-type stars. 

We attempted to constrain the source temperature by imple- 
menting a K-S test of the PI values of the observed events, as de- 
scribed above for arrival times. Model distributions of PI values 
were tested for an isothermal optically-thin plasma source, with an 
array of temperatures, T, at intervals of 0.25 dex in log T in the 
range 6.0 < logT (K) < 8.0 (Fig.|3ji)- Each source model was 
generated in XSPEC (Arnaud 1996) using an APEC model (Smith 
et al. 2001) with solar abundances and absorbing column density 
typical for the Pleiades, iVH= 2.0 x 10^" cm"^ (Stauffer 1984; 
Paresce 1984), convolved with the pn response matrix and ancil- 
liary response at the source position. 

While it was possible to rule out models with T < 10®'^ K 
at the 90 per cent confidence level, due to the small number of 
counts and the hard spectrum of the significant number of back- 
ground events we were unable to exclude high source temperatures, 
even up to W^° K (Fig.|3ji). 

For a plasma temperature of 10^" K the sensitivity ratio of 
pn to MOS is 3.6 in the 0.3^.5 keV energy band and the MOS 
count-to-flux conversion factor is 6.4 x lO^^'^ ergcni"^ s~^per 
(0.3^.5 keV) cs~^. Thus, the total MOS-equivalent effective ex- 
posure time for Roque 14 was 1 17.6 ks, its MOS-equivalent count- 
rate was 2.1 ± 0.5 X 10~^ c s~^, and its time-averaged X-ray lu- 
minosity was 3.3 ± 0.8 x 10^'' ergs"\ with Lx/Lboi « 10~''^°^. 
The assumption of a plasma temperature of 10®'^ K would give 
Lx«3.5±0.8 X 10^''ergs-\ 



4.4 Upper limits to X-ray emission from undetected brown 
dwarfs 

No significant X-ray emission was detected from any of the remain- 
ing four brown dwarfs* . At the NIR position of each brown dwarf, 
we counted events detected by each EPIC instrument in the 0.3- 

1 .4 keV energy range'' . A radius of 8 arcsec was used to minimize 
counts from nearby sources to Roque 1 1 and Teide 1. The expected 
number of background counts was estimated from the value at the 
source position in the reconstructed image generated by EMLDE- 
TECT in the 0.3-1.4 keV band. This enabled us to account for stray 
counts from nearby sources, and we note that the respective back- 
grounds for Teide 1 and Roque 1 1 were thus 25 and 36 per cent 
higher than the values at those positions in the ASMOOTH back- 
ground map. Upper limits to the source counts at 95 per cent con- 
fidence were calculated using the Bayesian method described by 
Kraft, Burrows & Nousek (1991). 

These upper limits were corrected for the enclosed energy 
fraction at the source off-axis angle (0.40-0.46 for a radius of 8 
arcsec; Ghizzardi 2002) and converted to MOS-equivalent count- 
rates by dividing by the value at the source position in the 
MOS-equivalent EPIC exposure map. The assumption of a 10^° 
K plasma with solar abundances and A^h= 2 x 10^" cm"'^ 
gives a MOS count-to-flux conversion factor of 7.2 x 10~*^ 
ergcm"^ s"*per (0.3-1.4 keV) c s~* and upper limits to the 0.3- 

4.5 keV X-ray luminosity in the range 1.9-3.4 x 10^^ erg s~*, with 



-10" 



(listed in 



upper limits to Lx/iboi in the range 10 ^■^''- 
Table|2}. 

To optimise our sensitivity to the detection of X-ray emission 
from these brown dwarfs and put the tightest possible constraint 
on their mean X-ray luminosity, we constructed a composite image 
of all the available EPIC data in 1 x 1 arcmin squares centred on 
the NIR positions of the four brown dwarfs (Fig. 0)). The total 
MOS-equivalent exposure time for the composite brown dwarf was 
439.3 ks. EMLDETECT detected the source ?» 25 arcsec west of 
Roque 11, but no source at the position of the composite brown 
dwarf. We calculated an upper limit to the mean X-ray luminosity 



* A source detected by a previous analysis 15 arcsec from Roque 9 only in 
the 0.8-1 .5 keV band is not recovered in this procedure, and was in any case 
considered more likely to be the chance alignment of a background source 
(Briggs & Pye 2003b). 

5 For a T = 10'^° K plasma with solar abundances and N-n= 2 X 10^° 
cm~ ^ the 0.3-1 .4 keV energy band enables detection at a given signafnoise 
ratio against the observed background spectrum for the lowest total of 0.3- 
4.5 keV source counts. The required number of source counts is only slowly 
increasing with the upper energy bound up to 4.5 keV and we do not find 
better constraint of the Roque 14 source variability or plasma temperature 
using this stricter energy cut. 
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Table 3. X-ray and Ho properties of ultracool field dwarfs. Columns show: (2) spectral type; (3) bolometric luminosity in 10^" ergs~^(as published for 
VB 8, VB 10 and LP 944-20; as calculated from Baraffe et al. 1998 for age > 10^ yr and T^ft = 2475 K for LHS 2065); (4) emission type (flaring, F, or 
non-flaring, N); (5) logarithmic ratio of X-ray and bolometric luminosities; (6) reference for X-ray emission; (7) equivalent width in Ha in A; (8) logarithmic 
ratio of Ha and bolometric luminosities (calculated as in Mohanty & Basri 2003); (9) reference for Ha emission; (10) ratio of X-ray and Ha luminosities. 



Name 


SpT 


^bol 


Em 


log(Lx/ibol) 


Rcfx 


EWhc 


log(iHc./ibol) 


RcfHa 


Lx/Lho, 


(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


(9) 


(10) 


VB8 


M7.0 


2.40 


F 




-2.8 


Schmitt et al. 1995 


20.4 


-3.58 


Tinney & Reid 1998 


6.0 








N 




-3.5 


Fleming el al. 1993 


6.8 


-4.06 


Mohanty & Basri 2003 


3.6 








N 




-4.1 


Fleming et al. 2003 


3.7 


-4.32 


Tinney & Reid 1998 


1.7 


VBIO 


M8.0 


L70 


F 


> 


-2.8 


Fleming et al. 2000 
















N 




-4.9 


Fleming et al. 2003 


5.6 


-4.32 


Mohanty & Basri 2003 


0.3 








N 


< 


-5.0 


Fleming et al. 2000 


4.1 


-4.46 


Martin 1999 


< 0.3 


LHS 2065 


M9.0 


L20 


F 




-2.5 


Schmitt & Liefke 2002 


261.0 


-2.84 


Mi>rtin& Ardila2001 


2.2 








N 


< 


-3.7 


Schmitt & Liefke 2002 


9.0 


-4.30 


Mohanty & Basri 2003 


<4.0 








N 


< 


-3.8 


Schmitt & Liefke 2002 


7.5 


-4.38 


Miirtin& Ardila2001 


< 3.8 


LP 944-20 


M9.0 


0.57 


F 




-3.7 


Rutledge et al. 2000 
















N 


< 


-5.7 


Rutledge et al. 2000 


1.0 


-5.26 


Mohanty & Basri 2003 


< 0.4 



of the four brown dwarfs of 1.1 x lO'^^ erg s~^ using the method 
described above. This corresponds to a mean Lx/Lboi of 10"^'*. 
Deeper observations or a composite analysis of a larger sample are 
required to determine if the mean activity level of Pleiades brown 
dwarfs Ues well below the saturated level. 

Considerably lower 2-a (~ 95 per cent confidence) upper lim- 
its, to Lx in the range 5.2-9.5 xlO^® ergs~^, and to Lx/iboi 



in the range 10 



-10" 



, have been reported for these ob- 



jects (excepting Roque 9) by Neuhauser et al. (1999; henceforth 
N99) using a number of observations of the Pleiades by the ROSAT 
PSPC. All five brown dwarfs were included in 7 separate exposures 
longer than 1.5 ks of 4 different PSPC pointings in the ROSAT pub- 
lic archive; the total exposure time of each field ranged from 22.5- 
39.9 ks. They were best-observed in the "Pleiades Center" field, at 
off-axis angles in the range 14-24 arcmin. The longest of the three 
exposures of this field, 22.4 ks from a total of 35.4 ks, was not used 
in N99, probably due to a faulty aspect solution (Hodgkin, Jame- 
son & Steele 1995). We have recalculated 95 per cent confidence 
upper limits by applying the Bayesian method as described above, 
using the broad-band (0. 1-2.4 keV) images, background maps and 
exposure maps of the 6 remaining archival exposures, and merg- 
ing data from exposures of the same field. We used a variety of 
extraction radii in each field and further combined data from differ- 
ent fields obtained with similar enclosed energy fraction (calculated 
using Boese 2000, equation 9), to find the strictest upper limit for 
each brown dwarf. We were also mindful of avoiding nearby X- 
ray sources detected by XMM-Newton. Optimum enclosed energy 
fractions ranged from 0.35-0.75. We converted PSPC count-rates 
to unabsorbed fluxes in the 0.1-2.4 keV band using a conversion 
factor of 1.0 x 10~^^ ergcm"^ s~'^per cs"'^, appropriate for a 
10^° K plasma with low absorption, as used in N99. Fluxes in the 
0.3^.5 keV energy band would be 5 per cent higher. A distance of 
135 pc was used to calculate luminosities (N99 used 125 pc). 

The tightest upper limits to the X-ray luminosity we could 
apply were 1.2-2.5 xlO^** ergs"^ for Roque 9, 14, 11, and 12, 
and 4.3 x 10^* erg s~^ for Teide 1, which may be contaminated by 
stray counts from the X-ray-bright K5 Pleiad HII 1348, ~ 50 arcsec 
away (Briggs & Pye 2003a). Upper limits to I/x/^boi were in the 
range 10~^^-10~^^. Hence, we conclude that the iJOSAJ obser- 
vations were not sufficient to detect Roque 14 at the level detected 
here by XMM-Newton, and the current XMM-Newton observation 
places the strictest upper limits thus far to the X-ray emission levels 
of Pleiades brown dwarfs of spectral type M7.5-8. 



5 DISCUSSION 

5.1 X-ray emission from ultracool dwarfs 

We have detected X-ray emission only from the M7-type Roque 14 
at I/x/iboi ~ 10^^"^, and placed an upper limit to the mean 
emission level from the four later-type Pleiades brown dwarfs at 
Lx/Lhoi < 10""^''. This is consistent both with the pattern of 
emission from active main-sequence M stars, which scales with 
bolometric luminosity, and with the idea that persistent X-ray emis- 
sion levels fall around spectral type M8 as magnetic field dissipates 
more easily in cooler, more neutral lower atmospheres. 

We can make some comparison of the character of the mag- 
netic activity of Roque 14 with that of other low-mass stars based 
on the relative levels of chromospheric and coronal emission. Flem- 
ing (1988) has reported a mean ratio of Lx/Lua ~ 6.7 for a sam- 
ple of the active, frequently-flaring, dMe stars and Reid, Hawley 

6 Mateo (1995) have reported a value of ~ 3 for MO-6 field stars 
outside of flares. Unfortunately, Lua and Lx are both typically 
variable and simultaneous measurements are scarce. Observed val- 
ues and calculated ratios of Lx and Lua for the Ha- and X-ray- 
detected dwarfs VB 8, VB 10, LHS 2065, and LP 944-20 are listed 
in Table |3] We stress that these Lx/Lna ratios are not calculated 
from simultaneous measurements of Lx and Lua , but from a sim- 
ple ranking of the observed values of each. Nevertheless, a strik- 
ing change in non-flaring emission appears to take place between 
the M7 star VB 8, in which Lx> Lhq and the M8 and M9 dwarfs 
VB 10 and LP 944-20, in which Lx < Ln^ (noted by Fleming et al. 
2003). This suggests that the efficiency of persistent coronal heat- 
ing decays more quickly in response to the increasingly neutral at- 
mosphere of ultracool dwarfs than the efficiency of chromospheric 
heating. Within the same sample, excluding LP 944-20, transition 
region heating, inferred from C IV emission, appears to remain as 
efficient as chromospheric heating (Hawley & Johns-Krull 2003). 

The ratio of the single measurements of Lx and Lua from 
Roque 14 is 4.0. If Roque 14's measured chromospheric and coro- 
nal emissions are interpreted as persistent this suggests that its 
magnetic activity is of similar character to that of VB 8 and ac- 
tive main-sequence M dwarfs. Although we cannot exclude that 
Roque 14's observed X-ray emission is solely the result of a flare 
with decay time-scale ~ 5 ks, like the observed high-level X-ray 
emission from M8 old-disc star VB 10 and M9 dwarfs LP 944- 
20, LHS 2065 and RXS Jl 15928.5-524717, the high Ha emission 
level of Roque 14 supports a persistently higher level of magnetic 
activity than on these cooler dwarfs. 
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Roque 14 has been suspected to be a near-equal-mass binary, 
but Lx/Lua > 2 even if the X-ray emission is interpreted as com- 
ing from two stars . An observation by the Hubble Space Telescope 
NICMOS camera does not support Roque 14's binarity, finding no 
companion of comparable brightness at separations > 0.1 arcsec, 
or 13.5 AU (Martin et al. 2000). While interaction between two 
close binary components could be influential in the X-ray produc- 
tion mechanism, the emission is at a similar level to that produced 
by magnetic activity on single active MS M dwarfs. 

The observed strong Ha emission from the M7.5 Roque 12 
hints at activity similar to that of Roque 14. In contrast, the mea- 
sured Ha emission levels of the M8 Pleiades brown dwarfs Teide 1 
and Roque 11 are very like that of the M8 VB 10, and we spec- 
ulate that the magnetic activity of these objects is already in the 
regime where Lx< Lua and their persistent X-ray emission levels 
are more than an order of magnitude lower than that of Roque 14. 
Deeper X-ray observations are required to test this prediction. 



dwarfs are the same, and that Roque 14 and LP 944-20 are repre- 
sentative of their respective ages. There is a large spread in the ob- 
served Lx/I/boi values of young brown dwarfs, a significant num- 
ber being undetected at upper limits < 10^^ (e.g. Imanishi et al. 
2001). Roque 14 has one of the highest Ha emission levels among 
Pleiades brown dwarf candidates (Zapatero Osorio et al. 1997b), 
while LP 944-20 has one of the lowest among M9 field dwarfs (Mo- 
hanty & Basri 2003). So, while we may expect the X-ray emission 
level of brown dwarfs to decrease as they age from « 125 to ~ 320 
Myr, as apparently observed by comparing Roque 14 and LP 944- 
20, the effect is probably exaggerated in choosing these two objects 
as representative. Larger, and necessarily deeper, surveys of the X- 
ray emission of brown dwarfs in the field and in the Pleiades (and 
other clusters) are required to make further progress in understand- 
ing the evolution of the magnetic activity of substellar and ultracool 
objects. 



5.2 The evolution of X-ray emission from brown dwarfs 

Very young brown dwarfs in star-forming regions are now routinely 
observed to emit X-rays at levels of Lx/Lhoi ~ 10~''-10^'' that 
arise from hot (T > 10^ K) plasma (e.g. Imanishi et al. 2001; 
Preibisch & Zinnecker 2002; Feigelson et al. 2002). A growing 
body of evidence indicates that young brown dwarfs experience 
accretion and outflows just like young low-mass stars (Jayaward- 
hana, Mohanty & Basri 2003). Their X-ray emission is very likely 
to be produced by similar means to T Tauri stars. This is prob- 
ably largely coronal emission as a result of magnetic activity, as 
in older low-mass stars, and as young substellar objects may have 
photospheres as warm as 2900 K (e.g. Baraffe et al. 1998), their at- 
mospheric conditions are ripe for the efficient coronal heating seen 
on active MS M5-6 stars. However, Ha emission may arise pre- 
dominantly from material accreting on to the young brown dwarf, 
rather than from a hot chromosphere. X-ray emission levels appear 
to rise with Ha emission levels up to EWua ~ 20 A, as would 
be expected for Ha emission predominantly from a magnetically- 
powered chromosphere, but there is no detection of X-rays from 
a young brown dwarf with EWna > 30 A (Tsuboi et al. 2003). 
This is consistent with the scenario emerging from X-ray surveys 
of T Tauri stars wherein samples of stars showing signs of strong 
accretion, such as high EWho, appear to show lower levels of X- 
ray emission than weakly-accreting objects of similar mass (e.g. 
Flaccomio, Micela & Sciortino 2003). 

X-ray emission has been previously detected from just two 
brown dwarfs older than 5 Myr. While any accretion is expected to 
have ceased by the age of the 12 Myr-old TWA 5B, this object has 
a much lower mass than Roque 14, 0.014—0.043 Mq (Neuhauser 
et al. 2000), and is unlikely to be a good model for its youthful 
X-ray activity. Conversely, LP 944-20 appears to have a similar 
mass to Roque 14, 0.056-0.064 Mq (Tinney 1998). It has cooled 
to a spectral type of M9, so its lower atmosphere is highly neutral, 
which is probably the key reason why its persistent X-ray emission 
level is ~ 100 times lower than those of many very young brown 
dwarfs and the detected level from Roque 14. However, we should 
be wary of assuming that the activity levels of all coeval brown 



^ As Liict/Lhol is calculated here from EWna^ as a ratio with the con- 
tinuum, it should be httle changed whether there is one source or are two 
similar sources. 



6 SUMMARY 

We have observed five candidate brown dwarfs in the Pleiades with 
XMM-Newton, detecting X-ray emission from the M7 Roque 14. 
The low number of counts and significant contribution of back- 
ground counts prevent meaningful constraint of the temperature 
and exclusion of transient emission. Assuming a plasma temper- 
ature of 10^" K and nominal absorbing column to the Pleiades 
of 2 X 10^" cm^^, the time-averaged X-ray luminosity of Lx 
= 3.3 ± 0.8 X 10^^ crgs^^, and its ratios with the bolometric 
(Lx/Lboi ~ 10"^"''') and Ha (Lx/Lua ~ 4.0) luminosities re- 
semble those of active main-sequence M dwarfs, as have been ob- 
served on the M7 old-disc star VB 8. 

We have placed the tightest upper limits thus far on the X-ray 
emission levels of the four later-type Pleiades brown dwarfs: Lx 
< 2.1-3.8 X 10^^ crgs-\andLx/-Lboi < lO^^^^'-lO^^'^^ 

The XMM-Newton data do not exclude that Roque 14's ob- 
served high level of X-ray emission is solely the result of a 
flare with decay time-scale ^ 5 ks like the observed high-level 
emission from M8-9 field dwarfs VB 10, LP 944-20, LHS 2065 
and RXS Jl 15928.5-524717, but the high Ha emission level of 
Roque 14 is supportive of a persistently higher level of magnetic 
activity than on these cooler dwarfs. However, the similarity of the 
low Ha emission levels of Teide 1 and Roque 1 1 to that of the 
M8 old-disc star VB 10 prompts us to speculate that the persis- 
tent X-ray emission levels of M8 brown dwarfs in the Pleiades may 
be low, with Lx/L^oi < 10"^, like those of VB 10 and the M9, 
^ 320 Myr-old brown dwarf LP 944-20. 

Deeper X-ray observations, less contaminated by background, 
are required to confirm a persistent high level of X-ray emission 
from Roque 14, to assess its coronal temperature, and to probe 
the typical X-ray emission levels of Pleiades brown dwarfs. Co- 
ordinated programmes to study X-ray, Ha and radio emission from 
ultracool and brown dwarfs of a range of ages are required to en- 
able significant progress toward understanding the mechanisms and 
evolution of magnetic activity on these low-mass, cool objects. 
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